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Abstract 

We investigate the Brans-Dicke (BD) theory with the potential as cosmological model to explain 
the present accelerating universe. In this work, we consider the BD field as a perfect fluid with 
the energy density and pressure in the Jordan frame. Introducing the power-law potential and 
the interaction with the cold dark matter, we obtain the phantom divide which is confirmed by 
the native and effective equation of state. Also we can describe the metric f{R) gravity with an 
appropriate potential, which shows a future crossing of phantom divide in viable f{R) gravity 
models when employing the native and effective equations of state. 
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I. INTRODUCTION 



Supernova (SUN la) observations has shown that our universe is accelerating 111. Also 

fi fi r 

cosmic microwave background radiation [2j, large scale structure |3|, and weak lensing [4| 
have indicated that the universe has been undergoing an accelerating phase since the recent 
past. Although there exist a number of models explaining an accelerating universe, the two 
promising candidates are the dark energy of cosmological constant in the general relativity jsl 
and a modified gravitational theory such as f{R) gravity js-S]. 

Recently, there was an extensive study of dark energy models based on the Brans-Dicke 
(BD) theory interacting with the cold dark matter (CDM) j^. However, the total equation 
of state Wcs = Ptot/Ptot was mainly used to measure the evolution of the universe. In order 
to see how the BD field describes the accelerating phase, we have to introduce both the 
native and effective equations of state because the BD field is non-minimally coupled to 



gravity 
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111 ]. Especially, we need the effective equation of state to take into account 



the universe evolution properly because there always exists an interaction between the BD 



fluid and the CDM 



121 ]. For the interacting holographic dark energy models, there is no 



phantorn phase when using the effective equation of state 13|] instead of the native equation 



of state 



14| . For the brane interacting holographic dark energy rnodels, the effective equation 



of state was used to account the evolution of the universe [15|. The effective equation of 
state could read off from the Bianchi identity which provides a non-standard conservation 
law. 

On the other hand, f{R) gravity models have been extensively employed to explain 
the present accelerating universe. The observational data might imply the crossing of the 
phantom divide Wde = — 1 in the near past [l^. In this case, the crossing of the phantom 



divide could be resolved in the viable f{R) gravity models [17l-ll9|. Especially, we would 
like to mention that a general approach to phantom divide in f{R) gravity was investigated 
in 18[, where the scalar-tensor version of f{R) gravity was used to see the phantom divide. 
However, it was shown that any singular f{R) gravity may be done non- sin gula r 20| . More 
recently, consistent, viable and non-singular f{R) gravity was suggested in 21]. 

Interestingly, it was shown that the viable four f{R) models generally exhibit the crossing 



of the phantom divide in the future evolution 



22]. 



A common feature to all analysis was performed by mapping the Starobinsky model 23] 
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to a scalar-tensor theory of gravity. It seems that the metric f{R) gravity is equivalent 
to the BD theory with curd = 0, while the Palatini f{R) gravity is equivalent to the BD 
theory with wbd = —3/2 [7|. Despite its mathematical equivalence, two theories may have 



shown physically non-equivalence 



24| : super-accelerating phase in the BD theory describes 



decelerating phase in /(-R) gravity. Also, it was pointed out that the mapping seems to be 
problematic because the scalar potential defined by U (^(R)) = R^ — f{R) with $ = dRf{R) 



25 



ml 



induces a singularity in the cosmological evolution 

Before we proceed, we wish to mention the difference between Einstein and Jordan 
frames jo]. We consider the frame in which non-relativistic matter (CDM, baryons) obey 
the standard continuity equation with ~ a^^- This is the Jordan frame as the physical 
frame in which physical quantities are compared to observations. It is sometimes useful 
to introduce the Einstein frame where a canonical scalar field is coupled to non-relativistic 
matter directly 28j. Even though one considers the same physics in both frames, using 
different time and length scales may offer the apparent difference between the observables 
in two frames. 

In this work we investigate how the present accelerating phase is realized in the scalar- 
tensor theory (BD cosmology). We consider the BD field as a perfect fluid with the energy 
density and pressure in the Jordan frame. Introducing the power-law potential (154|) and 
the interaction with the CDM, we confirmed the appearance of phantom divide by using 
the native and effective equation of state. Especially, inspired by the work of Ref. 22|, 
we study the cosmological implications of the /(-R) gravity using the BD theory with an 
appropriate BD potential ( l86l) . which indicates a future crossing of phantom divide in viable 
f{R) gravity models when employing the native and effective equations of state. This shows 
a close connection between BD theory and f{R) gravities for explaining future crossing of 
phantom divide. In the BD approach, we find a singularity in the past evolution of the 
universe. Hereafter, we consider the metric f{R) gravity only and thus, we mean f{R) 
gravity by the "metric /(-R)" gravity. 
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II. BD COSMOLOGY WITHOUT A POTENTIAL 



For cosmological purpose, we introduce the Brans-Dicke (BD) action with a matter in 
the Jordan frame 



S 



d'^x^—g 



$i? - WbD T + Cm 



(1) 



where $ is the BD scalar, cjbd is the parameter of BD theory, and Cm represents other 
matter which takes a perfect fluid form. The field equations for metric g^i, and BD scalar $ 



are 



1 

2' 
SttG 



2ujbd + 3 

where the energy-momentum tensor for the BD scalar is defined by 



(2) 
(3) 

(4) 



and the energy-momentum tensor for a perfect fluid takes the form 



Tja, = Pm.g^iv + (Pm, + Pm)u^Uy 



(5) 



Pm (Pm) denote the energy density (pressure) of the matter and is a four velocity vector 

with UaU°' = 1. 

Considering that our universe is homogeneous and isotropic, we work with the flat 
Friedmann-Robertson- Walker (FRW) spacetime 



dr + a\t) dr^ + r\de^ + sin 



(6) 



In this spacetime, the first Friedmann and BD scalar equations take the forms 
1$ 



6 v<i>; 3 $ ' 

MG{pm - 3pm) 



(7) 
(8) 



2a;BD + 3 

where H = d/a is the Hubble parameter and the overdot denotes the derivative with respect 
to time t. Here we note that the case of ojbd = —3/2 is not allowed unless a radiation- matter 
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FIG. 1: The equations of state W^-^ scalar versus its parameter wbd- ^ monoton- 

ically decreasing function of ojbd, while VFr-tj is a monotonically increasing function of oj^d- The 



bound of W-^Y) is given by —1/3 < Wnr^ < 1 because of wbd > — 3/2 and W 



BD 



1 as wbd — ^ CO- 



At wbd = 0, one finds that WgD ~ ^bd blows up. 



with pm = Pm/'i comes into the BD theory. Regarding the BD field as a perfect fluid, its 



energy and pressure are defined by kinetic terms as 
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Pbd 



1 



l^BD 



6i/ 



+ AH 



I 



(9) 
(10) 



If one does not specify the parameter cubd? one cannot determine the BD equation of state 
exactly. However, the Bianchi identity of V ^G^^ = implies that there exists an energy 
transfer between BD fluid and matter 



Pbd + 3i/(pBD + Pbd) = ^^f- 



(11) 



This continuity equation play a crucial role because it shows manifestly the energy transfer 
between pbd and pm and, thus, it defines the effective equation of state. 



On the other hand, we consider action with a minimally coupled scalar ip 

1 



28 



S 



LievrG 

in the Einstein frame. The field equations for metric ^^j^ and a scalar ip are 

G,. 



' fii/ 



= 0. 



(12) 

(13) 
(14) 



Its energy density and pressure for ip are the same given by 

= Y = P^^ (15) 

which describe a stiff matter with = 1. Here we obtain canonical forms for and 
p^, in comparison with non-canonical forms of pbd (El) and Pbd in the BD-frame. A 
continuity equation of + 3H {p^ + p^) = [unlike ffTTj) ] leads to the 7/^-scalar equation f|T^ 
exactly as 

^ + 3HtP = 0. (16) 
Importantly, the Bianchi identity leads to the two conservation laws separately 

V^T^t = ^ + Si/^/' = 0, V''T;^ = 0. (17) 
Hence, we need the EOS 

W^ = P^, (18) 

whereas we do not need to introduce the effective EOS W^^ like fl44l) arisen from fllip . 

In the absence of matter, the BD scalar plays a role of kinetic matter. This kinetic matter 
evolves as the conservation law is satisfied by itself 

PBD + SH{pBD + Pbd) = (19) 
whose equation of state (EOS) is defined by 

IVbd - (20) 
Pbd 

The solution to the first Friedmann and BD scalar equations is given by 



3("BD+l)±V'3(2c^BD+3) IT v/3(2^bD +3) 

a{t) = t 3(3^BD+4) ^ = t 3^BD+* . (21) 

Plugging the above into (120]) . one finds its EOS as 



w± 3(a;BD + 2) ± 2V3(2u;bd + 3) _ ^ 3 

W^^ = — o with c^BD > -t:- (22) 

owbd ^ 

In the limit of ubd 0, VFj^j-, — )■ 4/0, while Wj^^ |. Their behavior is shown in Fig. 
1. Here we choose Wbd = W^bd EOS for the BD kinetic-matter. The EOS bound is 
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given by —1/3 < Wbd < 1- If one requires the condition ( !TT|) together with = 0, the 
only allowable solution is the case saturating the lower bound 

^BD = -l^ ^BD = -I (23) 

which corresponds the solution to the conformal relativity: a{t) ~ t, $ ~ l/t^,|OBD ~ 
l/a^,PBD = — Pbd/3. This case gives a zero acceleration of a = 0. Consequently, the perfect 
fluid interpretation of the BD scalar is valid only for ubd = — 3/2 0. 

When the CDM is present, we have to solve the different equations. In the FRW space- 
time, equations take the forms 

^^^^(PB..^), (24) 
if = -4vrG(pBD+PBD + ^ + |^), (25) 
$ + 3i/$ = ^""^ (p^ - 3pm) , (26) 

^l^BD + -J 

where pm is the CDM density given by 

P. = § (27) 

with the present dark matter density pj^. It is convenient to use new variables as 

$' H' 

where ' denote the derivatie with respect to x. Also we define the density parameters 
Using the relations 

$ dxd^ I . , 

^ = -^^^ = H<^, (30) 
dt dx ^ ^ ^ ^ 

^ = + ^'- A^) , (31) 

energy density and pressure are given, respectively, by 

Pbd = [^BD^^ - Qf] , (32) 

Pbd = [wBD^^ + 4ip-2Xip + 2 {ip' + (33) 
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The Bianchi identity (fTT]) takes into account the energy transfer between BD field and CDM, 
while the CDM evolves according to its own conservation law 

Prn + 3H{p^+p^) = 0. (34) 

Eqs.dM]), (ES]) and can be written as 

1 = l^BD + ^^m, (35) 

, 3 AttG , , 

2 ^ -^Pbd, (36) 

- Av9 + 3<^ + ^ _ ri - ^bd), (37) 

where we used the pressureless condition of pm = for the CDM. Solving Eg. ( 1371) for and 
inserting it into pbd leads to 



Pbd 



u;bd^' - + - — ^-^(1 - ^bd) . (38) 



167rGL 2u;bd + 3 

Substituting this into Eq. fl36|) . we find 

X 3 1 
^=2 + 4 



c^BDV^' - + ^ ^ , J l-r^Bo) . (39) 



A further relation is found to be 

^^BD = 7(wbd¥^^ - 6(p). (40) 
o 



Eq. (1371) can be rewritten as 

Let us plug A and I^bd into Eqs. (14T|) and solve it numerically with the initial condition. 
On the other hand, we obtain the native EOS for the BD fiuid 

PBD ^BDy^'-2V^+2^(l-^BD) 

Pbd WbdV^^ - 



Considering Eq. f[TT|) as 

PBD + 3i7 (1 + W^^^) PBD = 0, (43) 
we obtain the effective EOS 

3pBD^ 

^ ^ ^ ^BDy'-^^+2^(l-^BD) ^ 

3 wbdV'^ - 
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class tp 17bd ^bd ^bd 

t \ 1 __5wBD+6 1 2(ajBD+l) 

<^BD+1 6{li;bd + 1)^ 3((^bd + 1) 5wbd+6 

(b) 5±^^^S 1 1^3+3 W^^ 

(C) ^"^tr"""'^ 1 ^BD ^BD 



TABLE I: List of critical points without BD potential. 
The initial value for is determined by 



^(0) = 3±V3(2n'^„.eD + 3) ^^^^ 

with 1]bd(0) = VL%^. Here +( — ) sign correspond to increasing (decreasing) ^bd at x = 0. 
For — signature, its evolution induces a singularity. This behavior could be expected from 
the critical points obtained by solving the equation of = 0. The result is summarized in 
the Table [H These critical points indicate asymptotic behaviors in the far future and far 
past. In order to test whether each critical point is or not stable, we need to observe the 
signature of dip' /dip. If the signature is negative (positiven), it may be stable for the far 
future evolution (far past evolution). The viable parameter range for class (a) is found by 
requiring the condition of < ^bd < 1 as 

3 4 6 
wbd < --, ~ 3 ^ ^BD < --• (46) 

However, if we demand the positive-definite energy density for the BD fluid pbd > and 



the negative-definite pressure pbd < [lO| , the relevant range is determined solely by 



3 

c^BD < --■ (47) 

In this case, the native and effective equations of state take the bounds 

W^BD, W^l > (48) 

which means that the BD fluid without potential does not explain the future phantom divide. 
For classes (b) and (c), these are nothing new because we have the condition 

UJBB >-l^ Wbb = W^l > (49) 
which corresponds to the absence of the CDM as is shown in Eq. (l22|) . 



Therefore, we note that the role of BD scalar without potential (equivalently, k-essence 
with non-canonical kinetic term only) as a source generating the accelerating universe is very 
restricted because it can at most describe "Wbd(W^bd) ~ —2/3 acceleration" in the presence 



of the CDM. In the presence of matters 



10|, the BD scalar $ appears to interpolate smoothly 



between the matter-dominated and accelerating eras by speeding up the expansion rate of 
the matter-dominated era like (a(t) ~ t^/s ^ ^a^^^ ^ 2{ubd + l)/(3a;BD + 4) > 2/3)), while 
slowing down that of accelerating phase derived by cosmological constant to some degree 
like (a(t) ~ e^* — (1 + ;^t)(2'^BD+i)/2-j_ Hence, we have to include an appropriate potential 
to obtain the phantom divide of VFbd(W^bd) — 



III. 



BD COSMOLOGY WITH A POTENTIAL 



The action for generalized BD theory is given by 11 1 

1 



S 



'S'n 



(50) 



where S^a = J d'^Xy/—gCm is the action for the other matter of the perfect fluid type and 
t/($) is a potential for the BD scalar $. The equation of motions can be obtained as 



(51) 



where T^^ is given by 



pBD 



8ttG 



^ I V^<1>V.<1> - ^g,. iV<^A + ^ { V^V.* - g,.V'^ - g,y87rGUi<^)} 



(52) 



and TJ^ is the energy- momentum tensor Equation for $ is changed to be 



2t/($) - 



SvrG 



(53) 



2a;BD + 3V'' ' J 2a;BD + 3' 

In the context of dark energy, it is possible to construct a single scalar model (scalar-tensor 
theory) based on the BD theory with constant f] = —dlnU{(j))/d(j). This depends on the 
choice of the potential U{^). 
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Setting $ = F[ip) = e the exponential potential Uoe takes the power-law form 



(54) 



with the constant a = ri/2Q. In the FRW spacetime (E]), three equations take the forms 



H = -47rG(pBD - -Pbd - ), 



f 2C7(*) - *^ 

2a;BD + 3 V 



These equations are consistent with Ref. 29| but are slightly different from Ref. 



SttG 



2uj- 



-XPm. - 3p Ti 



BD 



ing the BD field as a perfect fluid, its energy and pressure are changed as 

2 



3, 



30 
13| 



Pbd 



1 



167rG 
1 

167rG 



IJ-'BD 



+ 4i7- + 2- 



urn 



(55) 
(56) 

(57) 

. Regard- 

(58) 
(59) 



The Bianchi identity leads to the same relation as in Eq.ffTTl). 

In order to solve Eqs. (1551) . ( 1561) and ( l57j) . it is convenient to introduce new variables as 



(60) 



Using ( 130|) and ( 131]) . energy density and pressure are expressed in terms of ^ and -j/^, respec 
tively, by 



Pbd 
Pbd 



[wbd^^ + 4(p - 2A^ + 2 (v?' + (p2) - 2 



Then, Eqs. (l55l) . ( I56l) and (1571) can be written as 



^BD + — 1 
3 47rG 
2 



A 



i72 



Pbd, 



(2 - a) ^/^ 



BDJ 



2a;BD + 3 2cjbd + 3 

Solving Eq. ( l65l) for y?', and inserting it into j>bd, we obtain the pressure 



if2 



Pbd 



2 ^ 4(2 - a) , 

WBD'y^ - 2(/} + 



6 



2cjbd + 3 2cc;bd + 3 



'1 



BDj 



(61) 
(62) 

(63) 
(64) 

(65) 
(66) 
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Substituting this into Eq. (16^ . we arrive at 



, 3 1 

^ = 2 + 4 



c^BDV^' - + 5^ "I V^ + o ^ , o (l - ^bd) - . (67) 



2ubb + 3 2a;BD + 3 
Importantly, we note that Qbd can be written as 

^^BD = I (wbd<^' - + 2^) . (68) 



Plugging this into Eq. (165|) leads to 

9 2(2 -a) 3(1-1]bd) , / v 

2a;BD + 3 2a;BD + 3 

On the other hand, from the definition of tp, we obtain a newly differential equation 

^' = (ay, - + 2A) ip. (70) 

Now we have to solve two coupled equations (l69l) and (1701) for and t/^ numerically with 
initial conditions. 

Considering f l6T]) and fl66|) . we obtain the native EOS for BD field with potential 

PBD ^BPy^^ - 2^ + jgg^^ + ^(1 - I^BP) - 2^ 

This might not be suitable for representing the true equation of state for the BD scalar 
because of the non-conservation of this fiuid ( fTTj) . We remind the reader that ( |TT|) could be 
rewritten as 

PBD + 3i/ (1 + Wfj,) PBD = 0, (72) 

which implies an effective EOS 

^ ^ ^ ^BD^^ - + jg^V^ + ^(1 - ^bd) - 2^ ^ 
3 tJBDV^^ — 6v9 + 2-?/; 

A typical solution is given in Fig. [2] for ubd = 0, showing that Qbd + = 1- We observe 

that there exists a phantom divide (VFbd = W^bd — confirmed by effective EOS 

In deriving this numerical solution, it was necessary to impose the initial condition at 
ao = l(a; = 0) as an input of the current observation data. It is impossible to construct 
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X 



FIG. 2: Time evolution of the BD cosmology with potential U = Uq^"' with wbd = 0: 
Wbd (magenta) , VFg^(cyan), Oin(blue), Obd (green) . The initial condition at x = is imposed 
by a = 1.449, V(0) = 2.253, ip{0) = ^ - 0.75 = 0.001 (corresponding to JI^d = 0.75), and 
$(0) = 1.0. 



nearly past, present, and future acceleration phases without fixing the initial condition with 
the current observation data. Usually, we need one initial condition to solve the first order 
differential equation. Most of cases are needed to specify t = as the initial condition. 
However, since we do not know the origin of the dark energy clearly, we could not use t = 
as the initial condition to solve f l69p and f l70p for the accelerating phase. 

Since the effective gravitational constant G"^^ = G/^ varies with time, it should satisfy 
the observed limits defined as 



cff 



G 
G^ 



-H^ < 10""yr-^ 



(74) 
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It means that 



Ho^iO) < lO'^V"', (75) 
and, finally 

ip{0) < 0.0013. (76) 

Here we used the present Hubble parameter of Ho = 77km/s/Mpc = 2.5 x 10~^®s^^ = 
7.88 X 10'"yr-i. 

Considering the density parameter of BD field as 

^^^ = ^ = 6 ' ^^=3^' ^^^^ 

Qbd is regarded as dark energy density parameter and its current value will be determined 
by observation as I^bd- Hence two initial values of ip and ip are not independent, but they 
are related as 

UBDV^ - 6(/^(0) + 2^(0) - 6fi°D = 0, (78) 
which gives us 



^(0) = 3±v/9-2a;BD(^(0)-l]°J _ ^^^^ 



For ubd — 0, we can choose — sign so that 

^(0) = ^ - n^B^. (80) 



In order to find the asymptotic values for variable, we need to determine the critical points 
from Eqs.( l69|) and ( ITOl) . We list the critical points and corresponding physical variables in 
Table M and Table IHT] for lobd = 0. 

Now we wish to analyze which critical point is stable with time evolution. To determine 
the stability, we consider the perturbation around the critical points. If the coefficients 
of the perturbation is negative (positive), then it may be stable for future evolution (past 
evolution). Actually, classes (a), (b) and (c) are nothing new because these are exactly the 
same classes in the absence of the potential (see Table [T]). The class (d) is less interesting 
because both its native and effective equations of state do not provide the phantom divide. 
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class 




.1/1 

V 






T/f/ 
V^BD 


(aj 


1 


n 
U 


5wbd+6 


1 


2(tJBD + l) 




6(wbd+1)2 


Slti'BD+l) 


5(^BD+6 


(b) 


3+V6WBD+9 


U 


-1 

i 


^BD 


^BD 


(c) 


3- V6<^BD+9 





1 


W^D 


W^b"d 


(d) 


_3 


3(a+3wBD+3) 


7a+6a;BD+3 


_1 


2a 


a 


2a^ 


2a^ 


a 


7a+6a;BD+3 


(e) 


2(a-2) 


(2i^BD+3){6u;BD-(a-5)(a+l)} 


1 


2q:^— 9a— 6ci;bd+1 


2a^— 9a— 6cjbd+1 


a+2a;BD+l 


(q+2wbd+1)^ 


3(a+2wBD+l) 


3{a+2wBD+l) 



TABLE II: List of critical points with the power-law potential ([5 



class 






f^BD 




Wbd 


(a) 


1 





-1 


1 
3 


1 

3 


(b) 


N/A 










(c) 


-1 





1 


1 

3 


1 

3 


(d) 


_3 

a 


3{a+3) 
2a2 


7a+3 
2a^ 


_1 
a 


2a 
7a+3 


(e) 


2(a-2) 


3(a-5) 


1 


2a2-9a+l 


2a2-9a+l 


a+1 


(a+1) 


3(a+l) 


3(a+l) 



TABLE III: List of critical points in the presence of BD potential (j54p with liJbd = 0. 



The graph of equation of state Wbd as a function of a is shown in Fig. [3] with ubd = 0. 
The graph of effective equation of state Wby) as a function of a is given in Fig. H] with 
i^BD = 0. The class (b) is not available for ubd = 0, as was mentioned by the BD kinetic- 
matter using H^BD- 

The class (e) is a newly interesting case. Let us study the class (e) more. The solution 
to the phantom divide of W^^{a) = — 1 is given by 



«! = 1, 02 = 2, 

while the solution to the dust matter of W^^{a) = takes the forms 

9-V73 9 + V73 
as = ^ 0.114, 04 = ~ 4.386. 

The solutions to the radiation of ^^§0(0^) = | is given by 

as = 0, ae = 5. 



(81) 



(82) 



(83) 
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1 



^ 0.5 




-0.5- 



FIG. 3: The graphs of Wbd as a function of a for the case of ^jbd = 0. Class (a): red; (b): N/A; 
(c): green; (d): dark yellow ; (e): blue. 



The minimum value of W^{ol) 



/2 13 
- - - ^ -1.067, 



(84) 



appears at 



On 



-1 + ~ 1.449. 



(85) 



Finally, we mention the a-dependent evolutions of two equations of state VTbd and W^. 
As is shown in Fig. [5l for a < 1, there is no phantom divide, while for a > 1, there is 
phantom divide. For a > amm, there are two crossings of = — 1. For 1 < a < ttmin, 



there is one crossing of 



eff 
BD 



1 and H^BD approaches de Sitter value of —1 for a = 1, 2. 
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1 



0.5 




-0.5- 



FIG. 4: The graphs of as a function of a for the case of wbd = 0. Class (a): red (it is 

overlapped with green); (b): N/A; (c): green; (d): dark yellow ; (e): blue. 

IV. BD COSMOLOGY AS f{R) GRAVITY 



The crossing of the phantom divide could be understood in the viable f{R) gravity 



models 



19| . Recently, it was shown that the viable four f{R) models generally exhibit 



the crossing of the phantom divide in the future evolution when using the EOS wbe = 
Pde/Pde 22 1 . Hence, it is very important to see whether the future crossing of phantom 
divide is available for the BD cosmology with the corresponding potential. Hence, we analyze 
the BD cosmology with the potential 



(86) 
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FIG. 5: Time evolution of WBD(dotted curve) and VFg^Dl^olid curve) for a = 0.500 (blue), a = 1.000 
(green), a = 1.449 (red as in Fig. [2]), a = 2.000 (magenta), a = 2.200 (cyan). 

which could be obtained by considering the equivalence between the f{R) gravity and the 
scalar-tensor theory (BD theory with potential and wbd = 0) [3l|- Actually, this potential 
is an approximation to Hu-Sawicki and Starobinsky f{R) models. Here two parameters C 
and p are chosen 



0<C<1, 0<j9<l. 
Previous analysis could be applied to here when defining a to be 

where a is determined by 

$(1 _ <|,)p-i 



a($) = Cp 



l-C(l-$)? 



(87) 



(88) 



(89) 
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class V? V ^ ^BD ^BD ^BD 

\A) \J o 1 a 1 3 3 1-C(1-<I>)P 3 3 1-C{1-<S>)P 
(b) 1 -1 



(c) -10 1 

(d) 4 15 1 



1 _ 1 

3 3 

1 1 

3 3 

1 1 

3 3 



TABLE IV: List of critical points with the potential ()86p . Here Z is the solution of Z{p — 1) = 
In ,, -2 

C[e^(p-2)-p]- 

Then, we could use the same equations in Sec. IIIII Since a is a function of $, we need to 
solve the differential equation for $. Hence the full equations to be solved are given by 

^' = V - 3^ + ^ (2 - «($)) ij + ^ (1 - I^bd) + A^, (90) 

V'' = - V? + 2A) (91) 

$' = (92) 

Relevant variables take the forms 
3 1 



. 4(2 6 , ^ , 1 

(Jbd = ^{uJbbv' 'ef + 2;p), (94) 

•^BD = 5 ^ -771 > (95j 

'^"^ " 3 + c^BDV^^ - + 2^- ■ 

In order to check that our system is working properly, we first recover Tsujikawa et 



al's result [3l| by taking wbd = 9998. 5((5 = 0.01). We have recovered their result of fig. 
1 correctly, which is shown in Fig. [6l This figure shows that for < C < 1, the matter- 
dominated phase with Wcs — is followed by the de Sitter phase with Wcs — —1. In contrast 
to this, our equations of states Wbd and W^bd show that the BD field evolves from a stiff 
matter with lyBD(W^BD) = 1 in the far past to de Sitter phase with Wbd(W^bd) = ~1 
the far future. We note that the phantom divide appears in the near future. Also, we wish 



to point out that the initial condition used in 3l| do not provide a proper evolution. This 



means that the density parameter Qbd of BD field becomes negative, showing an unphysical 
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case. This explains why equations of sate diverge at some points. However, our modified 
initial condition gives a correct evolution for all relevant physical variables. Our initial 
condition is chosen by requiring the nonnegative density parameter of f^BDl^^min) > on 
whole evolution. In Figure O we have used Xmin = and ^BD{xmm) = 10"^. We note that 



the definition of energy density and pressure for the BD are different from those in [3l|, but 
the definition of density parameter ^bd is the same. Furthermore, they have analyzed the 
evolution only for the positive x > (future direction), as one can see from the Figure [61 
However, the general tendency of its whole evolution is almost the same as in the power-law 
potential f/ = f/o$". 

Furthermore, we have to mention that two equations of state Wud and H^bd defined by 
Eq.f lTT]) and (175]) are slightly different from the notation used in 31| 



^eff =-^-l§-2 (97) 

which was defined from the total conservation law 

Ptot + 3iJ(ptot + Ptot) = ptot + 3iJptot(l + Wcs) = 0, ptot = Pm + Prd (98) 

with pbd = '^'Pbd- Hence Wes represents the equation of state for the whole matters in the 
universe. 

The asymptotic critical points are listed in Table IIVI Actually, tow classes (b) and (c) 
are nothing new because these are nearly the same classes in the absence of the potential 
(see classes (a) and (c) in Table UTTl) . The class (d) is less interesting because it is similar to 
class (c) in Table UTTl) . We note that all of these belong to unstable critical points. On the 
other hand, the class (a) represents the asymptote of f{R) gravity models and its critical 
point indicates the par future behavior, showing that IVbd(W^bd) — — 1- 

In order to have the correspondence between f{R) gravity and the BD theory, we have 
to choose the BD parameter to be zero (wbd = 0). In this case, it is not easy to impose 
appropriate initial conditions of Q^n = 0.25 and Q^d = 0.75 by adjusting two potential 
parameters C and p. A time evolution is depicted in Fig. 13 As this figure is shown, there is 
a singularity at a past of x = —3.16 which refiects that the mapping is problematic. In this 
case, however, the equation of state Wbb approaches —1 (de Sitter spacetime) oscillatory as 
the universe evolves toward the far future. Also, the effective equation of state does 
show the nearly same behavior toward the far future. This confirm the presence of future 



crossing of phantom divide which appeared in the four viable f{R) models 22| clearly. 
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FIG. 6: The graphs of ^^(blue), Obd (green) , (magenta) , VFBD(brown), tt;cfr(red: definition 



used in Ref. j3l[|) as a function of x with wbd = 9998.5. The imposed condition is C = 0.7, Q = 0.01, 
p = 0.2, (p{0) = 6.47x10-^, = 2.85x10-3 (corresponding to Obd(O) = 10-^, $(o) = 5.0x10-^3^ 



V. DISCUSSIONS 



We have employed the BD cosmology to explain the accelerating universe and future 
crossing of phantom divide. In this work, we regard the BD field as a perfect fluid model. 
First, the role of BD scalar without potential (k-essence with non-canonical kinetic term 
only) as a source generating the accelerating universe is very restricted because it could 
describe 'Wbd{W^%) = -2/3 acceleration" in the presence of the CDM. 

Turning on the power-law potential ( l54l) . the BD cosmology could describe the accelerat- 
ing universe in the interaction with the CDM. In this case, we have used both the equation of 
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FIG. 7: A time evolution of BD fluid with potential ([86l) and ujbd = (/(i?)-gravity): 
VFbd (magenta) , VFg^(red), f]m(blue), Obd (green). The initial condition at x = is imposed 
as C = 0.82, p = 0.09, ip{0) = 2.35, V'(O) = 9.3 (corresponding to Q^^^ = 0.75), ^>(0) = 0.6. 
A singularity appears at x = —3.16 when evolving toward the far past, while there is a damped 
oscillatory evolution toward the far future. 

state Wbd and effective equation of state M^bd check whether the phantom phase appears. 
Exphcitly, the BD field acts as a radiation field in the far past, whereas it plays a role of 
phantom field in the far future. This is compared with the case without the potential where 
the BD field acts as a radiation field in the far past, while it plays a role of an accelerating 
matter in the far future. This shows that the presence of the potential is crucial for obtaining 
a phantom divide in the BD cosmology. 

Concerning the BD description of /(-R)-gravity, we have chosen the potential in Eq. (l86ll 
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inspired by Hu-Sawicki and Starobinsky f{R) models. For the case of ubd = 0, the evolution 
of Qrn. and Qbb are similar to the power-law potential fl5^ . but there exists a singularity 
at X = —3.16, which restricts evolving toward the far past of a; = — oo after imposing the 
initial condition at x = 0. However, we have found that the universe evolves toward the far 
future of X = oo nicely. Both native and effective equations of state converge to —1 (the de 
Sitter spacetime of class(a) in Table \W\ oscillatory, which indicates that the BD description 
is working for showing a future crossing of phantom divide appeared in the viable f{R) 



gravity models \22 \. 
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